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Abstract 

We investigate the possibility that inhomogeneous nucleosynthesis may even- 
tually be used to explain the abundances of 6 Li, 9 Be and B in population II stars. 
The present work differs from previous studies in that we have used a more ex- 
tensive reaction network. It is demonstrated that in the simplest scenario the 
abundances of the light elements with A < 7 constrain the separation of inho- 
mogeneities to sufficiently small scales that the model is indistinguishable from 
homogeneous nucleosynthesis and that the abundances of 6 Li, 9 Be and B are then 
below observations by several orders of magnitude. This conclusion does not de- 
pend on the 7 Li constraint. We also examine alternative scenarios which involve a 
post-nucleosynthesis reprocessing of the light elements to reproduce the observed 
abundances of Li and B, while allowing for a somewhat higher baryon density 
(still well below the cosmological critical density). Future B/H measurements 
may be able to exclude even this exotic scenario and further restrict primordial 
nucleosynthesis to approach the homogeneous model conclusions. 
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1. Introduction 



There has been considerable recent interest in the possibility that baryons may have 
been distributed inhomogeneously in the early universe. There are a number of mecha- 
nisms by which such inhomogeneities could be produced (c.f. Malaney & Mathews, 1993). 
Perhaps the most frequently considered has been a first-order QCD phase transition. It 
is quite possible that density inhomogeneities could be produced (Crawford and Schramm 
1982; Hogan 1983; Witten 1984; Iso, Kodama & Sato 1986; Fuller, Mathews & Alcock 
1988; Kurki-Suonio 1988; Kapusta and Olive 1988). These perturbations may have had a 
profound effect on the production of the light elements in the early universe (Applegate and 
Hogan 1985; Sale and Mathews 1986; Applegate, Hogan & Scherrer 1987; Alcock, Fuller 
& Mathews 1987; Malaney and Fowler 1988; Kurki-Suonio and Matzner 1989 and 1990; 
Terasawa and Sato 1989a,b,c; Kurki-Suonio, Matzner, Olive & Schramm 1990; Mathews, 
Meyer, Alcock & Fuller 1990). In this paper we examine the consequences of inhomoge- 
neous nucleosynthesis on the intermeditate stable isotopes 6 Li, 9 Be, 10 B and 11 B. 

In recent years measurements have been made in population II stars of elements 
that were once thought to have been produced in insignificant quantities in the standard 
homogeneous big bang. One of these, 9 Be (Rebolo et al. 1988; Ryan et al. 1992; Gilmore, 
Edvardsson, & Nissen 1991) has been proposed (Boyd, & Kajino 1989; Malaney, & Fowler 
1989) as a potential signature of baryon-inhomogeneous nucleosynthesis. However the 
strong observed correlation of Be/H with metallicity implies that it was made in the early 
Galaxy (Walker et al. 1993, Fields, Schramm & Truran 1993). To date, population II B/H 
data includes only three measurements (Duncan, Lambert, & Lemke 1992), which again 
show a correlation with metallicity. There has been one observation so far of 6 Li (Smith, 
Lambert, & Nissen 1992). 

These data appear to be best explained by galactic cosmic ray spallation (Steigman 
& Walker 1992; Walker et al. 1993; Prantzos, Casse, & Vangioni-Flam 1993; Olive and 
Schramm 1992; Steigman et al. 1993), however there remains the question of whether 
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6 Li, Be and B can be produced by primordial nucleosynthesis, and if so whether they can 
provide the much sought after litmus test of baryon inhomogeneities in the early universe. 
Four of us (Thomas et al. 1993, hereafter TSOF) used the largest network to date (in 
terms of reactions influencing light, A < 12, element abundances) to demonstrate that 
standard homogeneous big bang nucleosynthesis underproduces these nuclei by at least 2 
( 6 Li), 4 ( 9 Be) and 5 (B) orders of magnitude, when the abundances of the lighter elements 
and 7 Li are compared to the observations. In this paper we investigate inhomogeneous 
yields using an even further extension of the reaction network developed in TSOF. Our 
work primarily differs from previous similar studies (e.g. Terasawa & Sato, 1989c) in that 
we have considered a larger network of reactions to produce 6 Li, 9 Be, 10 B and 11 B, fully 
allowing for neutron-rich flows and multiple back-reactions. 

Nucleosynthesis in a homogeneous big bang requires the evolution of a set of equations 
representing the rates of the nuclear reactions in the network. The only input parameter 
is the baryon to photon ratio rj = ub/tIj or equivalently the density of baryons, (since n 7 
can be directly related to the microwave background temperature). Additional parameters 
are introduced when the effects of baryon inhomogeneities are taken into account. One of 
these parameters is the length scale, /, associated with the fluctuations. If this length scale 
is much greater than the neutron diffusion length then diffusion is unimportant and the 
yields can be obtained by simply averaging the yields from regions with different baryon 
densities (Wagoner, 1973; Yang et al. 1984). In the homogeneous model however, data 
on D and He restrict the baryon density to a small range about 7] 10 = 3 (7710 = 1]/I0~ 10 ) 
where 7 Li takes on its minimum value ( 7 Li/H ~ 10 -10 ) in agreement with observations 
(Spite & Spite 1982a,b,1986; Hobbs and Duncan 1987; Rebolo, Molaro & Beckman 1988). 
Any average of the rjio = 3 solution with that for another value of rj will increase the 7 Li 
abundance, and is likely to violate the upper limit (1.4 x 10~ 10 ). Thus one obtains strong 
constraints on the amplitude of such perturbations. 

In the other extreme (inhomogeneities much smaller than the neutron diffusion length) 
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baryon diffusion will eradicate the inhomogeneities before nucleosynthesis begins, and we 
return to the homogeneous case. It is the intermediate case that interests us here. When 
the inhomogeneity scale is of the same order as the neutron diffusion length the more rapid 
diffusion of neutrons (compared to protons) leads to an inhomogeneity in n/p in addition 
to the density inhomogeneity. The earliest studies of this scenario (Applegate, Hogan & 
Scherrer 1987; Alcock et al. 1987) assumed that neutrons diffused to a homogeneous density 
before nucleosynthesis began, and neglected all diffusion effects during nucleosynthesis. 
This simple model was able to satisfy constraints from 2 H and 4 He (and 3 He) with a 
baryon density equal to the critical value (Qb = 1) and a density contrast (~ 100) that 
seemed not implausible from the point of view of quark-hadron physics. This scenario was 
attractive because it did away with the need for non-baryonic dark matter. Unfortunately, 
it overproduced 7 Li. It was later suggested that the excess Li may be removed by diffusion 
of neutrons back into the high density region after nucleosynthesis (Malaney and Fowler 
1988) for values of / ~ 10 m. Though this turns out to be only partially true (Terasawa 
& Sato, 1989; Kurki-Suonio & Matzner, 1989; Mathews et al., 1990), the important lesson 
is that an accurate determination of abundances requires a calculation which takes careful 
account of the diffusion of neutrons and protons before, during and after nucleosynthesis. 

In recent years, more detailed diffusion calculations (Kurki-Suonio and Matzner 1989 
and 1990; Terasawa and Sato 1989a,b,c; Mathews, Meyer, Alcock & Fuller 1990; Kurki- 
Suonio, Matzner, Olive, & Schramm, 1990) have shown that not only could Li be affected 
but 4 He as well. It was found that nucleosynthesis with Qb = 1, no matter what the density 
contrast, overproduced both 4 He and 7 Li.^ Indeed these latest calculations all showed that 
for 0_b = 1, and when the distance scale of the inhomogeneities, / is greater than 30 m only 
the D abundance can be brought into agreement with observations. Though the standard 

t The problem with 4 He is particularly important since it does not allow for the possi- 
bility that consistency of all the light elements is achievable simply by the depletion of 7 Li 
in non-standard solar models 
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model constraints on r\ can be modified, the modification was shown to be rather limited 
(Kurki-Suonio, Matzner, Olive & Schramm 1990). 

In this paper we have used the diffusion code developed by Mathews et al. (1990). 
Initial density fluctuations are arranged in a lattice of spheres with separation / = 2r. Each 
sphere is described as a high density core and a low density outer shell. The core has density 

1 /3 

and radius of R and f v ' respectively, relative to the outer region. The sphere is divided 
up into concentric spherical zones, with a higher resolution near the boundary between 
high and low density regions. We expect the choice of spherical boundary conditions to 
maximize the potential effect. In all results presented here we have used 16 concentric 
zones. We have also run test cases with 8 and with 32 zones, indicating that 16 provide 
adequate accuracy, while remaining economical with computer time (Mathews et al. 1990). 
We also note that results of this code are consistent with those of Terasawa and Sato 
(1989a-c) and Kurki-Suonio and Matzner (1989). 

To alleviate some of the problems encountered in inhomogeneous models, mechanisms 
have been proposed to reprocess the nucleosynthesis products subsequent to the epoch 
of primordial nucleosynthesis. These mechanisms in particular reduce the abundance of 
7 Li, thus (potentially) allowing for higher fi B models. One such mechanism (Alcock et 
al., 1990; Jedamzik, Fuller & Mathews, 1993a) examines a fluid mechanical property of 
the electromagnetic plasma near the end of nucleosynthesis. The photon mean free path 
A 7 and the average physical (i.e. not comoving) size lh of the high density regions have 
different temperature dependences. At high temperatures (T > 20 keV), A 7 < lh, and 
so the EM plasma is confined over regions smaller than the baryon fluctuations, thus 
preserving these fluctuations. Below T = T m ~ 20 keV, however, A 7 > lh, and the EM 
plasma is not confined on the fluctuation scales. Protons in the high density regions flow 
out, hindered only by radiation (Thomson) drag. Alcock et al. (1990) argue that the 
dissipation of the fluctuations will homogenize the universe. They model this effect by 
running the inhomogeneous code to a given T m , then following the rest of the evolution 
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in the standard (homogeneous) code. They find that for the favored range of T m there 
is a significant reduction in the final abundances of Li, Be, and B over inhomogeneous 
production without dissipation, with 7 Li/H in particular reduced from ~ 10 -9 to ~ 10~ 10 . 
However, detailed calculations of Jedamzik et al. (1993) have shown that this mechanism 
is not as efficient as was previously believed. 

Gnedin and Ostriker (1993) have suggested another model of reprocessing, in which 
a baryon rich universe (Q° B ~ 0.15) overproduces 4 He and 7 Li, while underproducing 2 H 
and 3 He. They then posit that Jeans-mass black holes are formed at recombination. The 
black holes form accretion disks which emit a photon flux and reprocesses the ambient 
material; in particular, photodissociating the light elements and producing 2 H and 3 He by 
dissociating 4 He. The net effect could be to reproduce the observed levels of 2 H and 3 He, 
while still overproducing 7 Li by a factor of 10, and producing 4 He at a level of Y p ~ 0.250. 

We will discuss the implications of reprocessing on our conclusions. 

2. The calculation 

Based on the reaction network developed earlier (TSOF) and extended where neces- 
sary to allow for neutron-rich flows etc., we have evaluated the yields of the light elements 
from inhomogeneous primordial nucleosynthesis. The earlier reaction network contained 
180 reactions, see table 1 of TSOF. However, as pointed out in TSOF, the "flow" to the 
heavier elements (Be, B) lies largely along the neutron-rich side of the network, and thus 
mainly occurs in the low density, neutron rich zones of the inhomogeneous model. We 
felt it wise therefore to update the network further. The 84 additional reaction rates were 
estimated using the methods outlined in TSOF and Fowler and Hoyle (1964), and are 
shown in table 1. The full network is shown in figure 1. We have run a few sample cases 
without the extra reactions, and find no significant effect on the results. 

The diffusion code of Mathews et al. (1990) includes full multi-zoning and neutron 
back-diffusion. The diffusion coefficients used were those calculated by Banerjee & Chitre 
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(1991), and Kurki-Suonio et al. (1992). Results were obtained for a wide range of values 
for rj and r with the density contrast fixed at R = 10 2 , 10 3 , and fractional volumes (for 
the high density region) f v = 1/64, 1/8. We have also calculated results for R = 10 6 with 
f v = 1/64. The geometry assumed was spherical, with 16 concentric zones and a high 
density core. Due to the small uncertainty in the neutron mean life (± 2.1 sec) we fix this 
value at 889.1 s. (Particle Data Group, 1992). 

Baryon inhomogeneities have been best motivated by a possibility of a first-order QCD 
phase transition. Though the values of R and r can not yet be reliably predicted by QCD, 
some estimates can be made. For example, if chemical equilibrium is maintained, the value 
of R, which is very sensitive to the transition temperature T c ,is found to be between 7 < 
R < 100 for T > 100 MeV (Alcock et al., 1987; Kapusta & Olive, 1988). More generally, 
the value of R is determined by a combination of the enhanced thermodynamic solubility of 
baryon number in the high-temperature phase and the limited baryon number permeability 
of the moving phase boundary. Depending on the efficiency of baryon transport and the 
baryon penetrability of the phase boundary, R may be considerably larger (Witten, 1984; 
Fuller et al., 1988; Kurki-Suonio, 1988). The ultimate value at the time of nucleosynthesis, 
however, is expected to be less than 10 6 due to the effects of neutrino-induced heating 
and expansion of the fluctuations (Heckler & Hogan, 1993; Jedamzik et al., 1993a,b). The 
baryon number build-up at the boundary surface (where R is largest) contains only a small 
fraction of the total baryon number (Kurki-Suonio 1988). Thus, though we include values 
of R as large as 10 6 in our calculations, this should be viewed as an extreme upper limit. 

The value of r is also very sensitive to T c and the surface tension, a, of the phase 
interface (Fuller, Mathews & Alcock 1988; Kajantie, Karkkainen & Rummukainen 1990); 
r~2x 1 4 m( M ° v3 ) 3/2 (jf§v)~ 13/2 ■ For values of a 1/s ~ 70 MeV estimated by Fahri and 
Jaffe (1984) which agree with the effective field theory model estimates (Campbell, Ellis 
& Olive 1990), r < 0.4m for T c > 100 MeV. This is to be compared with preferred values 
of r ~ 30 m or the more recent estimates (Kurki-Suonio et al. 1992) of r ~ 100 m at 
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which reductions (though still insufficient) in the production of 4 He and 7 Li occur. It is 
important to note that the available estimates from QCD are all perfectly compatible with 
homogeneous nucleosynthesis. 

3. Results 

Results are shown in figures 2-8. Figure 2a shows the rj-r plane (where rj is the baryon 
to photon ratio, ns/n 7 and r is the radius of the spherical regions in cm, measured at 
100 MeV, after the phase transition) for R = 100 (results for f v = 1/8 and f v = 1/64 
are combined). The contours show observational limits on the abundances of the light 
elements (Walker et al. 1991 (WSSOK) and refs. therein): 

0.22 < Y p < 0.24 (1) 

2 H/H > 1.8 x 10" 5 (2) 

( 2 H+ 3 He)/H < 1.0 x 10" 4 (3) 

1.0 x 10" 10 < 7 Li/H < 1.4 x 10" 10 . (4) 

In addition, the dashed curve represents a He mass fraction Y p = 0.245 which is the 
most recently derived (preliminary) upper limit on Y p (Skillman et al. 1993). The region 
which satisfies all these constraints is hatched. Note that the only effect of increasing the 
maximum He abundance to 0.245 is to allow a slightly higher value of r. Figure 2b shows 
similar data for R = 1000, and 2c for R = 10 6 . Note that in figures 2 the hatched regions 
cover a similar area. 

For small r (< 100 cm), diffusion eliminates inhomogeneities before nucleosynthesis 
begins and the results are identical to those from a homogeneous calculation. As r increases 
from 100 cm, the He mass fraction rises rapidly above 0.24. Since all curves in figures 2 
are parallel to the r axis for r ^ 100 cm, we conclude that any inhomogeneous model 
that satisfies the limits on light element abundances will give the same abundances as 
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the homogeneous model, and the same limits on rj (2.8 < 7710 ^ 3.3, WSSOK). Since the 
hatched regions cover an almost identical area, this conclusion is independent of R and f v . 
We have also verified that in the large r limit (no diffusion) yields become independent of 
r. 

If we relax the upper limit on 7 Li, (say, because of some subsequent Li destruction) 
there is little change unless we also relax the upper limit on 4 He. The dashed curve in 
figures 2 represent a 4 He mass fraction of 0.245. In the case where Y p < 0.245 there are 
two allowed regions if we allow the primordial abundance of 7 Li to exceed 4 x 10~ 10 : (1) 
the previous limits are now 2.8 < i]io < 6 and r < 100 cm; (2) there is a region between 
the 2 H+ 3 He curve and the dashed curve at rj ~ 7, r ~ 10 4 . This solution however requires 
a rather finely tuned value of r in addition to the excess production of 7 Li, which would 
require the depletion of 7 Li by more than factor of 4 (we note that standard stellar models 
(Deliyannis, Demarque & Kawaler 1990) do not deplete Li significantly and non-standard 
stellar models which do deplete Li are highly constrained by the observation of 6 Li in HD 
84937 (Steigman et al. 1993)). Furthermore, since the code calculates abundances for a 
uniform lattice of spheres an accurate determination of yields in this case would require an 
averaging over a distribution of values for r. Given the narrowness of the allowed r values 
for the second solution it seems highly unlikely that realistic averaging would result in a 
solution satisfying all the light element abundances (Meyer et al. 1991) notwithstanding 
the problem with 7 Li. 

Abundances of 6 Li, 9 Be, 10 B and n B are shown in figures 3-6. With the exception of 
9 Be these are maximal abundances for all values of R, f v . Curves are given for 7/10 = 3.0, 
7.0, and 70.0. Yields of these elements are again independent of r for r < 100 cm, indicating 
that the yields are unchanged from those of homogeneous nucleosynthesis. For 7710 = 3 
this gives a 6 Li abundance (number density relative to H) of roughly 3 x 10~ 14 , a factor of 
100 lower than the recent measurement of Smith, Lambert, & Nissen (1992). Allowing for 
a higher 4 He abundance and abandoning the 7 Li constraints (that is, the r]io ~ 7, r ~ 10 4 
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solution mentioned earlier) increases the yield by (at most) a factor of 10. Of course if we 
now require the depletion of 7 Li, 6 Li will be severely depleted (Brown and Schramm 1988) 
and the discrepancy is amplified. 

The 9 Be abundance (figure 4) is shown as maximal abundances for R = 100, 1000, 
f v = 1/64, 1/8 (solid curve) and as abundances for R = 10 6 , f v = 1/64 (dashed curve). In 
this case, the effect of increasing R was greatest. For r < 100 cm, the 9 Be abundance is 
1 x 10~ 18 , four orders of magnitude below the observations (Rebolo et al. 1988; Ryan et 
al. 1992; Gilmore, Edvardsson, & Nissen 1991;). Allowing for 7710 ~ 5, r ~ 10 4 cm raises 
this almost to 10 14 (higher if we accept a density contrast of R = 10 6 ), however we regard 
this as an extremely unlikely situation. Note that even though 9 Be reaches a maximum 
at a few xl0~ 14 for 7710 = 70, the other light elements are irreconcilably off from their 
measured abundances, and this case is thus not viable. 

10 B (figure 5) and n B (figure 6) have abundances of ~ 10~ 19 and ~ 10~ 17 , respectively 
7 and 5 orders of magnitude below observations (Duncan, Lambert, & Lemke 1992). (The 
10 B abundance is always negligible compared to the 11 B abundance which is problematic 
if one wishes to show that the observed B is primordial; the two isotopes are observed with 
comparable abundances.) Using 7710 ~ 5 and r ~ 10 4 raises both abundance by less than 
two orders of magnitude. We emphasize that high abundances of 11 B are produced only 
for large 77, and at a value of r where the 9 Be abundance is low. Thus it is not possible 
to reconcile the large r model with the observed B to Be ratio, regardless of the problems 
with the light elements. 

We also show maximum (solid) and minimum (dashed) yields of 4 He and 7 Li in figures 
7 and 8. The possibility of a low 4 He abundance in inhomogeneous models was also 
investigated in Mathews, Schramm, & Meyer (1993). 
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4. Post-Processing 



We implemented the hydrodynamic-Thomson-drag dissipation effect by running the 
inhomogeneous code to a mixing temperature T m = 20 keV, the favored value given in 
Alcock et al. (1990). We then homogenized the results and continued to run down to the 
usual final temperature Tf = 10 7 K = 1.2 keV. While the effects of post-processing on the 
lightest element (D, 3 He, 7 Li) depend rather strongly on the input parameters, we did find 
consistent results for Be and B. The result is a reduction in the yield of 9 Be and 10 B, and 
an increase in n B. While it is conceivable that the right set of parameters might bring 
this model into agreement with the observations of the very light elements, the increase in 
B combined with a decrease in Be is difficult to reconcile with the observations (Duncan 
et al., 1992). Consequently, we regard this as an unlikely scenario. Similar conclusions 
have been reached in the recent parameter-free hydrodynamic calculations of Jedamzik & 
Fuller (1993). 



5. Conclusions 

With the accepted limits on the light element abundances ( 4 He, 2 H, 3 He, 7 Li) the 
length scale of inhomogeneities at the epoch of primordial nucleosynthesis is constrained 
to be r ^ 100 cm (at 100 MeV). With this constraint, the abundances of the light ele- 
ments, and of the additional elements 6 Li, 9 Be, 10 B and n B are largely indistinguishable 
from those of homogeneous nucleosynthesis. In particular, the abundances of LiBeB are 
lower (by several orders of magnitude) than the lowest of the abundances seen recently 
in population II halo stars. We conclude that these elements must be produced by some 
process other than primordial nucleosynthesis. 

If we push the limits on the light element abundances to the extreme, we find that 
while the abundances of LiBeB all increase, only 9 Be is raised significantly and still falls 
short of being able to explain any of the recent observations. 
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Table 1 

Reactions changed since TSOF 



Reaction Rate 



12 Be(/T) 12 B 34.31 ±0.03 

13 B(/T) 13 C 40.00 ±0.40 

14 B(/T) 14 C 43.05 ±3.21 

8 He(/T) 8 Li 5.68 ±0.09 

n Li(/r) n Be 79.7 ±0.9 

15 B(/T) 15 C 63 ±6 

17 C(/T) 17 N 3.43 ±0.29 

18 C(/T) 18 N 10.5i 2 ' 4 

18 N(/T) 18 1.10 ±0.06 

19 C(/T) 19 N 57.8 

19 N(/3") 19 3.3 ±1.7 

20 C(/T) 20 N 74.5 

20 N(/3") 20 O 6.9^; 4 

21 0(/3 ) 21 F 0.202 ±0.006 

22 0(/T) 22 F 0.308 ±0.67 

22 F(/T) 22 Ne 0.1639 ±0.0015 

23 F(/T) 23 Ne 0.31 ±0.06 

24 F(/T) 24 Ne 2.04 ±0.48 

18 0(n, 7 ) 19 21 ± 1 + 7.3 x 10 5 T~ 3/2 exp(-1.846/T 9 ) + 1.3 x 10 5 

n B(n, 7 ) 12 B 1.3 x 10 5 T 9 " 3/2 exp(-0.2112/T 9 ) ±4.0 x 10 5 T 9 " 3/2 exp(-4.53/T 9 ) 
±1.3 x 10 5 

12 B(n, 7 ) 13 B 1.3 x 10 6 T 9 " 3/2 exp(-1.64/T 9 ) ±1.3 x 10 5 

13 B(n, 7 ) 14 B 2.8 x 10 4 T 9 " 3/2 exp(-4.02/T 9 ) ± 1.3 x 10 5 

19 O(n, 7 ) 20 O 4.6 x 10 6 T 9 " 3/2 exp(-0.186/T 9 ) ± 4.6 x 10 6 T 9 " 3/2 exp(-1.74/T 9 ) 
±1.3 x 10 5 

21 F(n, 7 ) 22 F 1.8 x 10 6 T 9 " 3/2 exp(-4.18/T 9 ) ± 1.3 x 10 5 

n Be(n,7) 12 Be 1.3 x 10 5 

14 B(n, 7 ) 15 B 1.3 x 10 5 

16 C(n, 7 ) 17 C 1.3 x 10 5 

17 C(n, 7 ) 18 C 1.3 x 10 5 

18 C(n, 7 ) 19 C 1.3 x 10 5 

19 C(n, 7 ) 20 C 1.3 x 10 5 

17 N(n, 7 ) 18 N 1.3 x 10 5 

18 N(n, 7 ) 19 N 1.3 x 10 5 

19 N(n, 7 ) 20 N 1.3 x 10 5 

20 O(n, 7 ) 21 O 1.3 x 10 5 

21 0(n, 7 ) 22 1.3 x 10 5 

8 He(p,n) 8 Li 0.2874 x 10 11 T~ 2/3 exp(-6.4847/T 9 1/3 ) 

8 He(p, 7 ) 9 Li ditto 
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Table 1 (contd.) 



Reaction Rate 
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Table 1 (contd.) 



pa ct inn 


Rate 








19 N(a,p) 22 
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10 18 Tg 


~ 2 ^ 3 exp(- 
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20 N(p,n) 20 O 
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Figure Captions 



1 The nuclear reaction network used in all calculations. 

2a Limits on r and r\ due to the light element abundances, for R = 100. Curves show 
the most generous limits for f v = 1/8 and f v = 1/64, and represent the following 
abundances: 2 H/H = 1.8 x 10" 5 , ( 2 H+ 3 He)/H = 1.0 x 10" 4 , 7 Li/H = 1.4 x 10" 10 , 
Y p = 0.22, 0.24. The dashed curve is for Y p = 0.245. The hatched area shows the 
region allowed by the light element abundances. 

2b As figure 2, for R = 1000. 

2c As figure 2, but for R = 10 6 , f v = 1/64. 

3 6 Li abundance as a function of r for 7710 = 3, 7, 70 (7710 = 70 is included only for 
illustrative purposes, as the light element abundances can never all be fit in this case). 
The curves represent the most generous abundances for all values of R, f v . The 
hatched line shows the upper limit on r. 

4 As figure 3, with 9 Be abundances, except that the solid curves represent the maximum 
yields for all of R = 100, 1000, f v = 1/8, 1/64 and the dashed curves represent yields 
for R = 10 6 , f v = 1/64. 

5 As figure 3, with 10 B abundances. 

6 As figure 3, with 11 B abundances. 

7 As figure 3, with 4 He abundances. In addition, the dashed line shows the lowest yield 
for all .R, /„. 

8 As figure 7, with 7 Li abundances. 
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